Out-of-surface vortices in spherical shells 
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The interplay of topological defects with curvature is studied for out-of-surface magnetic vortices 
in thin spherical nanoshells. In the case of easy-surface Heisenberg magnet it is shown that the cur- 
vature of the underlying surface leads to a coupling between the localized out-of-surface component 
of the vortex with its delocalized in-surface structure, i.e. polarity- chirality coupling. 
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Understanding of the interplay between geometry and 
topology of condensed matter order is of fundamen- 
tal importance in several physical and biophysical con- 
texts, and this combination raises a number of unsolved 
questions. Exa mpl es are thin layers of superfluids and 
superconductors^^, nematic liquid crystal shells^, viral 
shells^, and cell membranes^. A considerable effort has 
been invested in understanding the role of the coupling 
between in-surface order and curvature of the underlying 
surface!^ The topological defects of 2D in-surface vec- 
tor fields are characterized by a winding number for the 
phase-like variable: vorticity g G Z (topological charge of 
7ri(S'^) homotopy group). On curved surfaces the Gaus- 
sian curvature leads to screening topological charges^. 
Vortices in curved superfluid films are a typical example 
of such kind of defects. 

Vortices in magnets belong to a more general type of 
topological defects. In addition to the vorticity, the mag- 
netic vortex is also characterized by the polarity p = ±1, 
which describes the vortex core magnetization. The topo- 
logical properties of magnetic vortices are characterized 
by the relative homotopy group 7T2{S'^^S^)^ and depend 
on both vorticity and polarity. Magnetic vortices were 
intensively studied during last decades for the sake of 
applications in nanomagnetism as high-density magnetic 
storage devices^^ and miniature sensor^^. Investigations 
of different aspects of magnetic vortex statics and dy- 
namics were mainly restricted to flat structures. In such 
nanomagnets, the vortex appears as a ground state in 
sub-micrometer sized magnets due to competition be- 
tween short-range exchange interaction and long-range 
dipole interaction^^ The ground state of smaller sam- 
ples is typically characterized by in-plane quasi-uniform 
magnetization. Contrary to in-surface, a quasi-uniform 
magnetization distribution in thin spherical shells is for- 
bidden for topological reasons; Instead, two oppositely 
disposed vortices are expected. 

In a flat nanomagnet, the vortex state is degenerated 
with respect to polarity. Hence, one can link the vortex 
polarity to the bit of information with possible spintron- 
ics application^^. One of the consequence of a more com- 



plicated topology of the magnetic vortex is a gyroscopical 
force that depends on both vorticity and polarity of the 
vortex. Therefore, the vortex polarity can be switched by 
exciting the gyroscopical motion. The switching thresh- 
olds for the two polarities are only equal in ideally flat 
structureJI^. Experiments on permalloy platelets have 
revealed a relatively large asymmetry in thresholdJ^^'^, 
which originates from the lack of the mirror-symmetry 
of rough thin-film structures'^. This indicates the inter- 
play between the vortex polarity and the curvature of the 
underlying surface. 

The influence of a curvature on magnetic proper- 
ties have been studied both experimentally and the- 
oretically for geometries of cylindei^'', torus'^, con^I^ 
and hemispherical cap structures'^. In this respect, we 
recently demonstrated experimentally by means of x- 
ray magnetic circular dichroism photoemission electron 
microscopy (XMCD-PEEM) the stability of magnetic 
vortex in thin permalloy (Ni8oFe2o) films on spherical 
particles'' . A precise theoretical description of pecu- 
liarities of vortices on spherical surfaces is not available 
in literature. Most theoretical studies are limited to 
skyrmion-like solutionJ^. 

In this paper, we study the structure of magnetic 
vortices on a thin spherical shell with an easy-surface 
anisotropy. Using anisotropic Heisenberg model, we find 
possible solutions of the vortex type. On the contrary to 
vortices in flat magnets, there is an interplay between the 
localized out-of-surface and the delocalized in-surface 
structure. In other words, the vortex core plays the role 
of a charge source for the vortex phase structure. 

The magnetic energy of a classical Heisenberg easy- 
surface ferromagnet has the following form: 



E = A dr 



-m • V'm 



(m • n) 



(1) 



with the exchange constant A, the anisotropy constant 
K > 0, magnetic length i = ^J AjK^ the surface nor- 
mal n, and the integration is over volume of the spher- 
ical shell. In the following, we use the local spherical 
reference frame for the unit magnetization vector m = 



2 



{rrir^m^^m^) = (cos 9, sin 6 cos ^, sin 6 sin Here, 
angular magnetic variables 6 = 6(r) and $ = <l>(r) 
describe the magnetization distribution with respect to 
the spherical coordinates (r, 'i^, x) of the radius-vector r. 
Hereinafter we consider a case of thin and high anisotropy 
shell: h <^ £ <^ where h is thickness of the shell and 
L is its inner radius. Therefore we assume that the mag- 
netization does not depend on the radial coordinate r. 
On the indicated conditions the total magnetic energy in 
terms of the local reference frame reads 
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E = AhJdxJ d^sin^l (S^e + cos^)^ 
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In this case, the static magnetization configuration the 
energy functional ([2| produces the Euler-Lagrange equa- 
tions 



V^e - sin 6 cos e 
sin^ 6 



(V<l>)^ - 1 + cot^ ^ 
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where S = cot cot — cos ^ and V-operators denotes 
the angular parts of the corresponding differential oper- 
ators in the spherical local basis. 

In the case of a high easy-surface anisotropy (^ ^ 0) 
the solution of Eqs. (|3| which minimize the energy ([2| 
reads O = 7r/2, <l> = const. This is a vortex solution 
where the magnetization is confined within the sphere 
surface except of two diametrically opposite point singu- 
larities - vortex cores. Such kind of "in-surface" vortices 
is well studied in different medicfe^. Here we demon- 
strate that taking into account the out-of-surface struc- 
ture of the vortex core (finite ^ > 0) essentially changes 
the vortex state properties in case of the curved surface 
as compared with planar magnets. In the following we 
consider only the azimuthally symmetric vortex solution 
= O('i^), ^ = ^{i^) by analogy with the planar vortices. 

The out-of-surface magnetization of the vortex core, 
so-called polarity, takes two values, p = ±1 (outward 
and inward). The magnetization distribution can be an- 
alyzed asymptotically near the vortex center {i} = 0), 
see Appendix [Aj The size of the vortex core is de- 
termined by the vortex out-of-surface magnetization, 
cosG ^ p{l - i9'^/2i9l), with i^c < 1. Moreover, the in- 
surface magnetization is described by the angular vari- 
able ^ $0 — psin$o'^^^/(4'i^c), where the constant 



will be determined later. Although the asymptotic limit 
of out-of-surface component is similar to that of a vor- 
tex in planar magnet^^, the phase ^ depends on the 
vortex polarity p, which is a distinct feature compared 
with the constant value in planar vortices. Since the out- 
of-surface magnetization has an exponentially localized 
structure, the following Ansatz function (similar to the 
vortices in planar magnets^^) can be used for description: 



cos 9 = Pie 2 
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where pi and p2 are polarities of the vortices at poles. 
Now, we consider the vortex in-surface magnetiza- 
tion. Accurate within the vortex core corrections, 
distribution can be described by the following equation: 



duu^ = -g{u)^\n^, g{u) 



(5) 



with u = lntan('i^/2) and mr{u) = cos 6. The function 
g{u) consists of two peaks localized near the vortex cores 
at Uc ~ lncot('i^c/2). 

In order to analyze the ^-distribution outside cores, 
we use the stepwise vortex shape model for the out-of- 
surface magnetization, rririu) ~ pi—pih{u-\-Uc)-\-p2h{u— 
Uc) with the Heaviside step function h{u). Using this ap- 
proach, g{u) becomes g{u) ^ ^^c[Pi^{u + Uc) — P2^{u — 
Uc)]. The consequence of such a model is that the local- 
ized out-of-surface structure plays the role of the charge 
density for the delocalized in-surface structure. The solu- 
tion of this model, which satisfies the Neumann boundary 
conditions du^{^oo) = 0, has the implicit form: 



^{u) = <^>o - -^^cPi sin ^{-Uc) [{u + Uc)^ - (^ - ^c)+] , 



pi sin<l>(-^Xc) =P2 sin^(iXc), 



(6) 



with ii+ = uh{u). The further analysis essentially de- 
pends on the relative orientations of vortices. 

For the case of same polarities {pi = P2 = p), the solu- 
tion <l> that minimizes the energy ^ takes the following 
explicit form outside vortex cores ydc <C 1): 



$(^) 
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where a is solution of the equation a = cos {ai^cUcTr/ 2). 
For details see Appendix [bJ Since 'dcUc <C 1, so a ^ 1. 
Accordingly, $ takes constant values inside the vortex 
cores, in particular, ^{i^ < i^c) ~ = i7r/2(l — 
pi^c^i^^c)- Energy of the vortex state with exception of 
the core energy Ec depends on the core size, E — Ec oc 
^luc. The dependence ^(d) is indicated in Fig. [l] by 
dashed lines. The approximate solution is in a good 
agreement with the numerical solution of Eq. ( |3b| ) , where 
the out-of-surface component was chosen according to 
Eq. Q. It should be emphasized that the phase of the 
vortex on a spherical surface gains a coordinate depen- 
dence given by Eq. ^ and has the maximum amplitude 
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FIG. 1. Possible vortex phases $ for the case of same po- 
larities pi — p2 — P (right inset). Soli d lines correspond 
to the exact numerical solution of Eq. (3b), where the out-of- 
surface component is chosen in form of Eq. J4| with 'dc — 0.05. 
The corresponding approximate solutions m are indicated by 
dashed lines. The corresponding magnetization distribution 
on the sphere surface is schematically shown in the left inset 
using arrows and stream lines. 




FIG. 3. Structure of the vortex state for the case of same po- 
larities (pi = p2 = 1) obtained by different methods. Line 1 — 
exact numerical solution of Eq. ^ with boundary conditions 
6(0) = e(7r) = and $'(0) = ^' {n) = and l/L = 0.05. 
Line 2 — (left) the Ansatz Q and (right) corresponding so- 
lution of Eq. (3b) where the function 0(1^) is determined by 
Eq. (|4| with an angular vortex core size 'dc — ^/ L — 0.05. 
This line coincides with line "p = 1" in Fig. ^ Line 3 and 
line 4 correspond to micromagnetic simulations of types (i) 
and (ii), respectively, (see text for details). Line 5 in the left 
graph shows the out-of-surface structure of the vortex core 
for the case of opposite polarities, when pi = —p2 = 1 and 

$ = TT. 
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FIG. 2. All possible vortex states of the spherical surface. The 
left column demonstrates the model out-of-surface magnetiza- 
tion distribution given by Eq. (|4| for all possible combinations 
of polarities. The right graph shows the corresponding dis- 
tributions of the phase $('i^). For the correspondence, we use 
the notation pip2, e.g. " + means pi = +1 and p2 = —1. 



in the center of each vortex, as opposed to the planar 
vortex. 

For the case of opposite polarities {pi = —p2 = p), 
the energy reaches its minimum for the trivial solutions 
^ = 7T for p = 1 and $ = for p = — 1 (as earlier, we 
consider the case i^c <C 1), see Appendix [B] for details. 
Such a solution can be considered as a three dimensional 
generalization of well-known onion state in narrow nanor- 
ings. We refer to this solution as the pumpkin state. The 
energy of the pumpkin state, E — Ec (x —"da is lower 



than for the vortex state. Namely, the energy gain is 
/S.E (X i^c- Nevertheless, it should be emphasized that 
these two states are separated by a high energy barrier 
related to the polarity switching of one of the vortices. 
Therefore, we suggest that both states can be realized 
experimentally even at room temperature. All possible 
vortex-like states described above are presented in Fig.[2j 

In this paragraph, a comparison between the obtained 
analytical results and the exact numerical solutions as 
well as the micromagnetic simulations is given. The nu- 
merical solution of Eq. d3| with the boundary conditions 
9(0) e(7r) and ^\{)) ^' {n) leads to a vortex 
structure for the case of same polarities pi = P2 = 1- The 
obtained out-of-surface structure of the vortex core cos 6 
is quite close to the model solution (|4| when dc = i/L 
(left graph of Fig. [3|. The solution for the vortex phase 
for same boundary conditions corresponds to two- 
fold degenerated solutions of the type given by Eq. ^ 
with opposite chiralities. One of such solutions is plot- 
ted in the right graph (Fig. |3|. The exact solution 
has slightly larger amplitude of turning compared with 
the model solution. This originates from a larger effec- 
tive core size i^c of the exact solution compared with the 
model given in Eq. (|4| (left graph of Fig. [3|. 

In order to verify our results, we performed two 
types of micromagnetic simulations by using the OOMMF 
code^^. A thin spherical shell was simulated considering 

(i) local magnetic interaction in the form of Eq. ([T]) and 

(ii) exchange and non-local magnetostatic interaction. In 
both cases, the material parameters^ were chosen in the 
way to provide the same ratio of characteristic magnetic 
length and the sphere radius £/L = 0.05. Physically, 
these two types of simulations are equivalent in case of 
vanishing thickness, when the magnetostatic interaction 
can be reduced to the easy-surface anisotropy. 

The simulations of type (i) confirm the analytical re- 
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suits with a high accuracy. In the case of same polarities, 
the system relaxes to the vortex state with the additional 
turning described by Eq. ^ (line 3 in Fig.|3|. For oppo- 
site polarities, the system relaxes to the state with ^ = tt 
or <l> = 0, as described above. 

According to the simulations of type (ii), in case 
Pi = P2i the magnetostatic interaction attenuates (but 
does not suppress) the phenomenon of the vortex phase 
turning (line 4 in Fig. |3| due to an increased energy of 
volume magnetostatic charges. In case of pi = — p2, 
the shell of mentioned size relaxes to the vortex state 
with ^ ^ zb7r/2 instead of the pumpkin state, which 
is preferred for spheres of smaller sizd^. To obtain 
an approximate criterion of the separation between vor- 
tex and pumpkin states, the difference of energies of 
the pumpkin and vortex states are estimated as follow- 
ing: AE = AEms + ^^ex^ where AEms ^ Lh^ is en- 
ergy increase due to volume magnetostatic charges and 
AEex ^ —i^ci'^h is the corresponding exchange energy 
decrease. Thereby, the pumpkin state is energetically 
preferable when L'^h < The detailed study of ground 
states of soft-magnetic spherical shells goes beyond the 
scope of this paper and it is subject of a prospective work. 

In conclusions, we predict novel features of magnetic 
vortex in thin spherical shell in comparison with the well 
known vortex in planar easy-plane magnet. We show 
that the vortex on a spherical surface gains a coordinate 
and polarity- dependent turning of its phase. An interplay 
between topological properties of the vortex, namely, its 
polarity, and the curvature of the underlying surface re- 
sults breaks the degeneration of the phase-like variable 
$ with respect to the rotation by any constant angle ^o- 
This degeneration is known as well as for tti -vortices in 
different media and for 7r2-vortices in flat magnets. It is 
instructive to note that the angle $o in magnetic nan- 
odisks determines the vortex chiralit}E3. Thus one can 
speak about polarity -chirality coupling. 



Appendix A: Asymptotic of functions S{i9) and 
in neighborhood of the vortex origin 

We consider here an asymptotic solution of Eqs. ([3| in 
neighborhood of point = 0. The functions 6(i?) and 
^{'d) can be presented in form of the Tailor series 



N 



e ^ 7rh{-p) + + ^ an^", (Ala) 

n=2 

N 



(Alb) 



where p is the vortex polar ity an d h{x) is the Heaviside 
function. The expansion ( |Ala ) satisfies the necessary 
boundary conditions: 6(0) = for p = +1 and 9(0) = tt 
for p = — 1. To obtain the asymptotic expansion accu- 
rate within terms of order (D{i9^) we can restrict ourselves 
with N = u where m = 2 is the order of Eqs. ([3|. 




Pl=-P2=^ 




FIG. 4. Function g(u) from Eq. ^ for = 0.05. 



Then we substitute the series (Al) into (|3|, expand the 
obtained equations by the small quantity ^ and equate 
the series coefficients of the same order terms until the 
order v. The obtained system results a relation of expan- 
sion coefficients in ( Al ). For the case v = 2 the described 
procedure results 



a2 =bi = 0, 

sin <I>o 



p- 



(A2) 



4^c 



Thereby we obtain the following asymptotic expansion 
accurate within terms of order 0{i^'^) 
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Appendix B: Vortex phase solution for the 

stepwise vortex shape model 

In the following, we focus on the solution of the Eq. ([5|. 
Using the stepwise vortex shape model mr{u)^ pi — 
Pih{u + Uc) + P2h{u — Uc)^ one can rewrite ^ in the 
form: 

duu^ = -9{u) sin<I>, 

TT (Bl) 

g{u) ^ -de [piK^^ + ^c) - P2^(^ - Uc)] , 

where 5{u) is the Dirac delta function. This approxima- 
tion agrees well with the real peaked form of ^-function, 
see Fig. [4j The general solution of (Bl) takes the form: 



Pi sm^{-Uc){u + Uc) + 
-P2 sin$(^Xc)(^ - ^c)+ 



(B2) 



with Uj^ = uh{u) and h{u) being the Heaviside step 
function. Using the Neumann boundary conditions 
9n^(±oo) = 0, one can find that = and 

Pi sin ^{—Uc) = P2 sin ^{uc). (B3) 

Substituting ([B3| into ([B2| results the Eq. ^ 



^{u) = ^0- TT^cPi sin ^{-Uc) [{u + Uc)^ 



{U-Uc) + ]. 

(B4) 
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Thereby for the interval —Uc<u<Uc one can write 
^{u) = au -\- where constants a and b can be found 
from the system 



Pi sm{—auc -\-b) = p2 sm{aUc + 6), 
a — ——VcPi sm[—auc + 6). 



The further analysis essentially depends on the relative 
polarities of vortices. For the case of equal polarities 
Pi = P2 = P the system (B5) has the following solutions 



where k G Z . Tak ing into account that 'dcUc ^ 1 one 
obtains from ( B6b ) only the trivial solution a = and b = 
kir. Substituting the obtained solutions into Hamiltonian 
([2| results that the trivial solution corresponds to the 
energy maximums and the solutions (B6a) minimize the 
energy. Therefore for the case of equal polarities one 
obtains the solution 

For t he c ase of opposite polarities pi = —p2 = p the 
system (B5) results 



a= — (2A: + 1), cos 6 = (-1)^(2A: + 1) — , (B7a) 

2Uc UcVc 



2p 



(B7b) 



b = ±|, a = T^i9cPcos{aUc), 
a = k7r, sin6= (-1)^+1-^, 



(B6a) 
(B6b) 



Due to the condition i^cUc <C 1 the system (B7) results in 
the trivial solutions a = and 6 = /ctt, and consequently 
<l> = or $ = TT. Analysis of the energy functional shows 
that for p = 1 the solution ^ = tt and for p = —1 the 
solution <l> = correspond to the energy minimum. 
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